Ion selectivity and current saturation in inward-rectifier K + channels 
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We investigated the features of the inward-rectifier K channel Kirl.l (ROMK) that underlie the saturation of cur- 
rents through these channels as a function of permeant ion concentration. We compared values of maximal cur- 
rents and apparent Km for three permeant ions: K + , Rb + , and NH 4 + . Compared with K + (i max = 4.6 pA and Km = 10 mM 
at —100 mV), Rb + had a lower permeability, a lower i max (1.8 pA), and a higher Km (26 mM). For NH 4 + , the 
permeability was reduced more with smaller changes in i max (3.7 pA) and Km (16 mM). We assessed the role of a 
site near the outer mouth of channel in the saturation process. This site could be occupied by either permeant ions 
or low-affinity blocking ions such as Na + , Li + , Mg 2+ , and Ca 2+ with similar voltage dependence (apparent valence, 
0.15-0.20). It prefers Mg 2+ over Ca 2+ and has a monovalent cation selectivity, based on the ability to displace Mg 2+ , 
of K + > Li + ~ Na + > Rb + ~ NH 4 + . Conversely, in the presence of Mg 2+ , the Km for K + conductance was substantially 
increased. The ability of Mg 2+ to block the channels was reduced when four negatively charged amino acids in the 
extracellular domain of the channel were mutated to neutral residues. The apparent K, n for K + conduction was 
unchanged by these mutations under control conditions but became sensitive to the presence of external negative 
charges when residual divalent cations were chelated with EDTA. The results suggest that a binding site in the 
outer mouth of the pore controls current saturation. Permeability is more affected by interactions with other sites 
within the selectivity filter. Most features of permeation (and block) could be simulated by a five-state kinetic model 
of ion movement through the channel. 



INTRODUCTION 

In addition to K + itself, inward-rectifier K + channels con- 
duct other ions including Rb + , Tl + , and NH 4 + . In Kirl.l 
(ROMK), Rb + behaves similarly to K + but has a lower 
conductance and permeability. NH 4 + is the only other 
physiologically important cation besides K + that is con- 
ducted by these channels. The selectivity properties of 
NH 4 + are interesting because the relative permeability 
PnhVPr measured under bi-ionic conditions is 0.1, 
whereas the relative conductance gNnVgic is approxi- 
mately one (Chepilko et al., 1995; Choe et al., 2000). 
This suggests that NH 4 + ions can pass through the chan- 
nels easily but have weaker binding interactions within 
the pore. K + could then readily displace NH 4 + from the 
pore, accounting for the low NH 4 + permeability under 
bi-ionic conditions. 

Inward-rectifier K channels have a high affinity for 
K + ions as assessed by the concentration dependence 
of conduction. In both Kirl.l and Kir2.1, currents sat- 
urate at low concentrations, with apparent K m values 
on the order of 10 mM (Lu and MacKinnon, 1994a; 
D'Avanzo et al., 2005). In the present study, we inves- 
tigated the relationship between Kirl.l cation affinity 
and permeability. In general, Kirl.l had lower affini- 
ties for less permeable cations. However, the relation- 
ship between conductance, permeability, and K m is 
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complex and appears to involve interactions with dif- 
ferent parts of the channel. Permeability depends on 
binding affinities within the selectivity filter, whereas 
K m is strongly influenced by sites outside the filter. 

MATERIALS AND METHODS 

Expression of Kir1 .1 in Xenopus laevis oocytes 

Oocytes were harvested from Xenopus according to the guidelines 
of, and with the approval of, the Institutional Animal Care and 
Use Committee of Weill Cornell Medical College. The animals 
were anesthetized through immersion in 1 liter of tap water 
containing 1.9 g L _1 tricaine methanesulphonate and HEPES 
(adjusted to pH 7.4) for 5-10 min. Once the animals were anes- 
thetized, a small incision was made in the abdomen and part of 
the ovary was removed. The oocytes were then dissociated 
through incubation in OR2 solution (in mM: 82.5 NaCl, 2.5 KC1, 

1 MgCl 2 , 1 Na 2 HP0 4 , and 5 HEPES, pH 7.4) supplemented with 

2 mg mP 1 collagenase type II (Worthington) and 2 mg mf 1 hyal- 
uronidase type II (Sigma-Aldrich) for 1 h. 

pSport plasmids containing Kirl.lb (ROMK2) were linearized 
with NotI restriction enzyme, and cRNAs were transcribed with 
T7 RNA polymerase using mMESSAGE m MACHINE T7 kit 
(Invitrogen). cRNA pellets were dissolved in nuclease-free water 
and stored at — 70°C before use. The oocytes were injected with 
10 ng RNA and incubated overnight in L-15 solution supplemented 
with HEPES, pH 7.4, 63 mg L -1 penicillin, and 145 mg L -1 
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streptomycin at 18°C. All chemicals were from Sigma-Aldrich 
unless otherwise noted. 

Patch clamp 

Before use, the vitelline membranes of the oocytes were mechani- 
cally removed in a hypertonic solution containing 200 mM su- 
crose. Patch-clamp pipettes were prepared from hematocrit 
capillary glass (VWR International) using a vertical puller (David 
Kopf Instruments) . They had resistances of 2-8 Mfl when filled 
with 110 mM KC1. For measurements in inside-out patches, pi- 
pette and bath solutions contained Cl~ salts of K + , Rb + , and NH 4 + , 
as indicated, plus 5 mM HEPES buffered to pH 7.4. In addition, 
bath solutions contained 0.5 mM EDTA. Currents from cell-at- 
tached and excised inside-out patches were recorded with a 
patch-clamp amplifier (EPC-7; HEKA), digitized with an interface 
(Digidata 1332A; Axon Instruments). Data were filtered at 1 kHz 
and analyzed with pCLAMP9 software (Axon Instruments). 

Data analysis 

Currents were measured from selected recording intervals normally 
containing a single active channel using all-points histograms. Histo- 
grams were fit with double Gaussian functions, and the current 
amplitude was defined as the distance between the two peaks. 

Single-channel currents (i) were analyzed according to the 
equation: 

i = i_/(l + K m /[c + ]), (1) 

where [C + ] is the concentration of the permeant ion. 

The effects of external blocking ions were analyzed according 
to the equation: 

i(B) = i(0)/(l + [B]/^), (2) 



with 

Kj(V) = Kj(0) • exp(z5FV/RT), ( 3 ) 

where [B] is the blocker concentration, z is the blocker valence, 
V is the transmembrane voltage, and 8 is the effective valence or 
fraction of the electric field at the blocking site. Kj(0) and 8 were 
estimated from linear regression analysis of a linearized form 
of Eq. 2: 

ln(i(0)/i(B) - 1) = ln([B]/Kj(0)) - zSFV/RT. 



Kinetic modeling 

The kinetic model shown in Fig. 10 was evaluated using Matlab 
7.8.0 as described previously (Edvinsson et al., 2011). Rate con- 
stants were calculated based on the values set for binding ener- 
gies, energy barriers, and electrical distances. The values for the 
fractional electrical distances 8i, &2, and 8 3 were 0.07, 0.113, and 
0.233, respectively, as used in previous studies (Kutluay et al., 
2005) . The model was fitted to the appropriate experimental data 
using the Isqcuruefit function in Matlab. 

Online supplemental material 

Online supplemental material shows the effects of external Na + 
on K + currents in the presence of the divalent cation chelator 
EDTA (Fig. SI). Figs. S2 and S3 illustrate the sensitivities of the 
measurable parameters of conductance, apparent Km, and bi-ionic 
reversal potentials simulated by these models to changes in the 
binding energies to the different sites. Table SI lists the kinetic 
parameters of the permeation model based on fits to experimental 
data. Figs. S1-S3 and Table SI are available at http://www.jgp 
.org/cgi/content/full/jgp.201110727/DCl. 
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Figure 1 . Conduction through 
Kirl.l channels in inside-out 
patches with symmetrical K + con- 
centrations. (A) Inward currents 
with —100 mV across the patch 
at different [K + ]. (B) i-V relation- 
ships for different [K + ]. Lines 
through the data points have no 
theoretical meaning. (C) Single- 
channel currents at —100 mV as 
a function of ion concentration. 
The line represents the best fit 
of Eq. 1 to the data, with i max = 
4.6 ± 0.1 pA and ^ = 10.2 ± 
1.2 mM. (D) i max and K m , mea- 
sured as in C, as a function 
of voltage. 
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RESULTS 

Conductance-concentration relationships 

We first confirmed the relationship between K + concen- 
tration and inward currents in Kir 1.1 channels. Chan- 
nels were expressed in Xenopus oocytes, and excised 
inside-out patches were formed with identical K + con- 
centrations, ranging from 5 to 200 mM, on both sides 
of the membrane. K + was increased by the addition of 
KC1 to the basic medium without cation substitution. 
Typical inward current traces for a membrane potential 
of —100 mV are shown in Fig. 1 A, and the corre- 
sponding i-V relationships are plotted in Fig. 1 B. At 
low K + concentrations, the open probability decreased 
sharply. Decreased open times, together with decreased 
currents, determined the lower limit of [K + ] at which 
currents could be resolved. The i-V relationships were 
approximately linear over the range of —20 to —100 mV 
at all concentrations. Plotting the current as a func- 
tion of ion concentration at —100 mV indicated a 
hyperbolic relationship with a maximal inward cur- 
rent (imax) of — 4.6 pA and an apparent K m of 10.2 mM 
(Fig. 1 C). The K m value decreased mildly with depo- 
larization of the membrane voltage, suggesting differ- 
ent voltage dependencies of rates of entering and exiting 



the channel (Fig. 1 D) . Outward currents were less reli- 
ably measured as a result of faster kinetics and in- 
creased noise levels. 

We repeated these experiments using Rb + as the con- 
ducted ion. Fig. 2 A shows inward currents at — 100 mV 
in inside-out patches with various Rb + concentrations 
from 10 to 200 mM on both sides of the membrane. 
Fig. 2 B shows i-V relationships for these concentra- 
tions. Rb + currents were lower than those for K + , con- 
firming previous findings (Chepilko et al., 1995; Choe 
et al., 2000) . At - 1 00 mV, i max was 1 .8 pA, whereas the ap- 
parent Km was 26 mM, higher than that for K + (Fig. 2 C) . 
The i-V relationship for Rb + was less linear than that for 
K + ; inward currents increased more rapidly with hyper- 
polarization. In contrast, the voltage dependence of the 
apparent K, n for Rb + conduction was less than that for 
K + (Fig. 2 D). The lower affinity for conduction corre- 
lates with the reduced permeability for this ion relative 
to K + , estimated from bi-ionic reversal potentials 
(Chepilko et al., 1995; Choe et al., 2000). 

We also examined the concentration-conductance 
relationships for NH 4 + . Handling of this ion is un- 
usual in that it has a permeability ratio Pnm/Pk of only 
M).l, determined from reversal potential measure- 
ments, but its conductance is comparable to that of 
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Figure 2. Conduction through Kirl.l channels in inside-out patches with symmetrical Rb + concentrations. (A) Inward currents with 
— 100 mV across the patch at different [Rb + ] . (B) i-V relationships for different [Rb + ] . Lines through the data points have no theoretical 
meaning. (C) Single-channel currents at —100 mVas a function of ion concentration. The line represents the best fit of Eq. 1 to the data, 
with i max = 1.8 ± 0.1 pA and Km = 26 ± 4 mM. (D) i max and measured as in C, as a function of voltage. 
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K + (Chepilko et al., 1995; Choe et al., 2000). We hypoth- 
esized that the low permeability might be explained 
if NH 4 + interacted less strongly (i.e., were bound less 
tightly) than K + as it traversed the channel, predicting a 
higher apparent Km for conduction. Current traces for 
NH 4 + are shown in Fig. 3 A, and i-V relationships are 
shown in Fig. 3 B. As found previously (Choe et al., 
2000) , the channel kinetics change when NH 4 + is con- 
ducted, with much briefer openings and lower open 
probability. With low NH 4 + concentrations, the open 
times were further decreased, limiting the lower end 
of the usable concentration range to 25 mM. As shown 
in Fig. 3 C, the maximal inward current for NH 4 + was 
3.7 pA at —100 mV, whereas the apparent Km value was 
16 mM, higher than but not strikingly different from that 
for K + . The voltage dependence of K, n was small, at least 
between —60 and —150 mV (Fig. 3 D). However, outward 
NH 4 + currents were virtually independent of NH 4 + con- 
centration over this range (Fig. 3 B), implying that the K, n 
for conduction in this direction was well below those for 
inward currents. However, as with K + , these currents were 
more difficult to measure accurately. 

To see if the saturation of NH 4 + currents is affected by 
the presence of K + on the opposite side of the mem- 
brane, we investigated inside-out patches in which NH 4 + 
was varied only in the pipette (extracellular) solution 
with constant 110 mM K + in the bath. Because the rever- 
sal potential changes under these conditions, we ana- 
lyzed inward NH 4 + slope conductances, rather than 
currents, over the concentration range of 50-200 mM 



(Fig. 4). These conductances (60-80 pS) were larger 
than those measured under symmetrical conditions for 
K + or NH 4 + . There was little change in conductance over 
this concentration range, again indicating high-affinity 
interactions with the channel. 

Occupancy of an outer binding site 

One possible site of interaction of permeant ions with 
the channel that could contribute to saturation of cur- 
rents is the outer mouth of the pore, just external to the 
selectivity filter. This could contain an obligatory bind- 
ing site for ions that controls the rate of movement into 
the channel. Analysis of a multi-ion pore with simplified 
association-dissociation kinetics showed that the pres- 
ence of such a site could lead to Michaelis-Menten-type 
kinetics of permeation (Nelson, 2011). Similar results 
are presented below for a more complex kinetic scheme. 

Because such a site might not be as strongly ion selec- 
tive as those of the filter itself, we examined whether the 
smaller alkali metal cations Na + and Li + could reduce K + 
conductance through the channel by competing for oc- 
cupancy. Fig. 5 (A and B) shows traces and i-V relation- 
ships in inside-out patches with 110 mM K + on the 
cytoplasmic side and either 1 1 mM K + alone or 1 1 mM 
K + plus 99 mM Na + on the extracellular side of the 
membrane. The effects of Na + were modest but clear. 
Outward currents were virtually unaffected, whereas 
inward currents were reduced by 50% or more. 

We analyzed these results using a simple Woodhull- 
type model (see Eqs. 2 and 3 in Materials and methods) 
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Figure 3. Conduction through 
Kirl.l channels in inside- 
out patches with symmet- 
rical NH 4 + concentrations. 
(A) Inward currents with 
— 100 mV across the patch 
at different [NH 4 + ] . (B) i-V 
relationships for different 
[NH 4 + ]. Lines through the 
data points have no theo- 
retical meaning. (C) Single- 
channel currents at — 1 00 mV 
as a function of ion concen- 
tration. The line represents 
the best fit of Eq. 1 to the 
data, with i max = 3.7 ± 0.1 pA 
and K„ = 16 ± 3 mM. (D) i max 
and Kn,, measured as in C, as 
a function of voltage. 
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that assumes a blocking site at a fixed point within the 
transmembrane electric field. The results of this analy- 
sis are shown in Fig. 5 D, in which the equation for frac- 
tional block is linearized with respect to voltage. The 
apparent Ki(0) obtained by extrapolating the linear 
regression line to V = 0 was 272 mM, whereas the slope 
of the line gave an effective valence of the blocking re- 
action of 0.15. That is, the voltage dependence of block 
could be explained by Na + blocking at a site that senses 
15% of the transmembrane electric field. In this inter- 
pretation, the effects of Na + appear as a reduction in 
single-channel current rather than a decreased open 
time because the blocking kinetics are too fast to resolve. 

Very similar results were obtained with Li + (Fig. 5, B 
and C). In this case, the estimated 1^(0) was 203 mM and 
the effective valence was 0. 1 7. The voltage dependence of 
the reduction suggests that Na + and Li + also bind weakly 
but specifically within the permeation pathway. 

Divalent cations interact with an outer 
binding/blocking site 

We next asked whether divalent cations could block the 
channel in a similar manner. Fig. 6 (A and B) shows 
traces and i-V relationships with 11 mM K + (and 99 mM 
Na + ) in the pipette in the absence of divalents and in the 
presence of Ca 2+ or Mg 2+ . Both Ca 2+ and Mg 2+ produced a 



voltage-dependent decrease in current, but the Mg 2+ - 
dependent block was stronger. Fig. 6 C shows the analy- 
sis of the voltage dependence; the K;(0) for Ca 2+ was 
about threefold higher than that for Mg 2+ , whereas the 
effective valences were similar. These values of 0.17 to 
0.18 are close to those measured for Na + and Li + , de- 
spite the fact that these blockers carry twice the positive 
charge. The voltage dependence of the block as well as 
the selectivity between these two divalent cations for 
block provides further evidence for a specific interac- 
tion with the pore. 

We extended these results by examining block by 
3 mM Mg 2+ at different extracellular K + concentrations. 
The effects of Mg 2+ are similar in each case, but the mag- 
nitude of current reduction increases as K + decreases. 
This is shown quantitatively in Fig. 7 A. The K ; (0) val- 
ues, interpolated with 1 1 mM K + and extrapolated with 
55 and 110 mM, increase with K + concentration, whereas 
the slopes of the plots are similar. The Ki(0) values are 
plotted as a function of K + in Fig. 7 B. The approxi- 
mately linear relationship is consistent with a competitive 
interaction between K + and Mg 2+ , with an estimated K; 
for Mg 2+ in the absence of K + of 2-5 mM and a dissocia- 
tion constant of 15 mM for K + binding to the Mg 2+ - 
blocking site. This is similar to the apparent for 
conduction (11 mM) in the absence of blocking ions, 
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suggesting that this site may well contribute to con- 
ductance saturation. 

We assessed the selectivity of this site by comparing 
the ability of different monovalent cations to relieve 
Mg 2+ block. This was assessed as the ratio of currents at 
-100 mV with and without 3 mM Mg 2+ and is plotted in 
Fig. 7 C. For Rb + and NH 4 + , external K + was completely 
replaced by these ions. For Na + and Li + , 99 mM of the 
test ion was added to a solution containing 1 1 mM K + . 
These results indicate that the site is somewhat selective 
for K + ; the other ions tested had similar abilities to dis- 
place Mg 2+ , with Na + and Li + being slightly more effec- 
tive than Rb + or NH 4 + . 

If this outer binding site determines saturation of 
currents, competition between K + and Mg 2+ in the 
outer part of the permeation path predicts that Mg 2+ 
should shift the apparent affinity of the channel for 
K + . To test this, we repeated the measurements of 
concentration-conductance relationships for K + in the 
presence of 3 mM Mg 2+ (Fig. 7 D). Under these condi- 
tions, the apparent K + affinity decreased markedly, 
with apparent Km values increasing from 11 to 96 mM. 
This is consistent with the idea that saturation of K + 



conductance involves K + occupancy of the outer bind- 
ing/blocking site, at least when Mg 2+ is present in the 
external solution. 

Role of fixed negative charges in saturation and block 
Affinities for ions at the outer mouth of the pore should 
be affected by the presence of fixed negative charges in 
the extracellular domain of the channel protein (D'Avanzo 
et al., 2005; Chang et al., 2010). To test this, we mutated 
several aspartate and glutamate residues in the outer 
mouth of the pore (E92, D97, E104, and E132) to eliminate 
their negative charges. The largest effect was observed 
with alterations at the El 04 position; this side chain is 
situated just above the outer aspect of the selectivity filter 
(Fig. 8 A). The E104S mutant significantly reduced the 
efficacy of Mg 2+ block without changing the voltage de- 
pendence (Fig. 8 B) . We also studied a quadruple mutant 
in which all four negatively charged amino acids in the 
extracellular domain were neutralized. Again, the affinity 
for Mg 2+ was decreased modestly but significandy, whereas 
the voltage dependence remained constant (Fig. 8 C). 

The elimination of negative charges did not have much 
effect on the apparent affinity for K + (Fig. 8 D) . This result 
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did not support the idea that the K m reflects binding 
to the outer mouth. However, we considered the possibil- 
ity that the charges might be screened, particularly by 
trace metal cations in the solutions. Indeed, inclusion of 
0.5 mM EDTA to chelate divalent cations significantly 
decreased the apparent Km to values too low to measure 
accurately (Fig. 8 D). The chelator did not change the Km 
in mutant channels lacking the negative charges in the 
outer mouth, consistent with the idea that the charges 
do affect affinity but are partially screened by residual 
divalent cations under most conditions. EDTA also 
modestly increased the block of the channels by ex- 
tracellular Na + (Fig. SI). 

In some Kir channels, fixed negative charges in the 
transmembrane cavity enhance the affinity for block 
by Mg 2+ and polyamines (Lu and MacKinnon, 1994b; 
Taglialatela et al., 1995; Yang et al, 1995). To see if such 
charges on the other side of the selectivity filter could 
also influence current saturation, we used the Kirl.lb 
mutant N152D, in which an aspartate is inserted in the 
second membrane-spanning segment at the position cor- 
responding to D172 of IRK1. Indeed, there was a signifi- 
cant decrease in the apparent K m for K + conductance, 



from 10.2 ± 1.2 (Fig. 1) to 4.7 ± 1.4 mM (Fig. 9). The effect 
is in the same direction as, but much smaller than, the in- 
crease in affinity for block by Mg 2+ or spermine block from 
the cytoplasmic side of the pore (Lu and MacKinnon, 
1994b; Taglialatela et al., 1995; Yang et al., 1995) . 

Kinetic modeling 

To examine if the saturation properties of the Kir 1.1 
channel could be accounted for with standard models 
of K + channel permeation, we used a five-state scheme 
used previously to describe KcsA channels (Kuduay et al., 
2005) and Kir4.1 channels (Edvinsson et al., 2011). The 
model can be formulated in terms of rate constants be- 
tween the states of the channel as described in Fig. 10 A. 
This scheme is based on known properties of the selectiv- 
ity filter of K + channels. Because Kir 1.1 and other inward 
rectifiers contain a cytoplasmic domain that may provide 
an additional resistance to ion flow (Choe et al., 2000) , 
the step between the cavity and the internal solution 
would include movement through this part of the pore. 

The model could fit data for i-V relationships with 
symmetric K + or NH 4 + concentrations, but many dif- 
ferent parameter sets gave equally good fits. We then 
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in the external solution. (D) Analysis of block by external Mg 2+ and Ca 2+ according to Eqs. 2 and 3. The lines represent best fits to the 
data, with Ki Mg = 13 mM and z8 Mg = 0.18, and Kc a = 31 mM and z8 Ca = 0.17. 
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attempted to simultaneously fit data from symmetrical 
K + , symmetrical NH 4 + , and bi-ionic K + versus NH 4 + con- 
ditions. We assumed the same kinetic scheme with the 
same number of sites and barriers with the same place- 
ment within the electric field for the two ions. Only the 
rates of transitions among the states were allowed to 
vary. A parameter set that gave satisfactory results is 
illustrated in Fig. 10 (B-E) . The pseudo-energy profiles 
for NH 4 + and K + along the pore for this set are shown in 
Fig. 10 F. These values are plotted as "rate-constant rep- 
resentation" (RCR) units (Andersen, 1999), arbitrarily 
assuming a frequency factor for the Eyring transition- 
state reaction rates of 6 x 10 10 s~ . The actual rates used 
are given in Table SI. 

These parameters were able to reproduce the basic 
properties of K + conduction as well as several of the 
remarkable features of the NH 4 + data, including the 
near independence of outward currents on ion con- 
centration, the high conductance to NH 4 + despite the 
low bi-ionic permeability, and the increased inward 
NH 4 + conductance with K + on the inside of the mem- 
brane. One prediction of the model that is different 
from the experimental observations is the slight de- 
cline in currents at high NH 4 + concentrations. The 
main features of the parameter set for NH 4 + are de- 
creased depths of energy wells for NH 4 + compared 
with K + , particularly at sites SI and S2. Although other 
parameter sets could also successfully fit the data, 
all of them involved similar changes at these sites. 



Thus, this simplified model of ion permeation can 
largely account for movements of K + and NH 4 + through 
the channel. 

We also performed a sensitivity analysis to examine 
the influence of model parameters on the predicted 
observables of apparent K m , conductance, and rever- 
sal potentials. When the binding sites were altered by 
simultaneously changing the energies of the well and 
adjacent barriers, the apparent Km values were most 
sensitive to increases in SO; increases of 1.3 RCR in- 
creased K, n by nearly 10-fold (Fig. S2 A) . Changes in SI 
had smaller effects, whereas Km values were almost inde- 
pendent of S2-binding energies. This may reflect the 
nature of the kinetic scheme, as the model is con- 
strained such that the S1-S4 sites' selectivity filter is 
occupied by two ions at all times. In contrast, inward 
conductance through the channel was most sensitive 
to changes in SI (Fig. S2 B) . 

Under bi-ionic conditions, with the parameters for 
the internal ion fixed (to those estimated for K + ) and 
those for the external ion varied, reversal potentials 
were inversely related to the binding energy at each of 
the sites. When all sites were raised or lowered simulta- 
neously, an increase of 1 RCR produced a reversal poten- 
tial of —60 mV (Fig. S3 A). These changes also reduced 
the open-channel conductance by ^50% (Fig. S3 B). 

The decreased energy wells for NH 4 + account for the 
decreased Pnhu/Pk ratio because under bi-ionic con- 
ditions, the sites will be preferentially occupied by K + . 
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The concomitant high inward NH 4 + conductance is a 
distributed property of the various energy barriers and 
is not so simply accounted for by any one parameter. 
We have not explored in depth the basis for other prop- 
erties, such as the higher inward NH 4 + conductance in 
the presence of internal K + or the low apparent K m for 
outward NH 4 + currents. 

DISCUSSION 

The main conclusions of this study are: (a) Currents 
through Kir 1.1 channels carried by NH 4 + saturate with a 
low apparent K m only modestly higher than that for K + , 
despite a reduced NH 4 + permeability, (b) Saturation of 
inward currents partly reflects occupancy of an outer 
binding site that is selective for K + but also accepts 



other monovalent and divalent cations, (c) The unusual 
permeation properties of the NH 4 + ion can be accounted 
for by conventional kinetic models of ion transport 
through Kir 1.1. These topics are discussed in more 
detail below. 

Saturation of K + currents 

We found that currents through Kirl.l channels satu- 
rated at low K + concentrations, with apparent Km values 
of ^10 mM. This is in general agreement with previous 
results on Kirl.l (Lu and MacKinnon, 1994a) and 
Kir2.1 (Lopatin and Nichols, 1996; D'Avanzo et al., 
2005; Chang et al., 2010), although the analyses of the 
concentration dependence varied. At very high K + con- 
centrations (>300 mM), conductances through Kirl.l 
channels can actually decline from maximal values 
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Figure 8. Removal of external 
fixed negative charges reduces af- 
finity for Mg 2+ block and K + per- 
meation. (A) Homology model of 
the outer mouth of Kirl.lb show- 
ing the positions of four negatively 
charged amino acid side chains, 
E92, D97, E104, and E132, rela- 
tive to the conduction pathway. 
(B) Analysis of block by 3 mM Mg 2+ 
of wild-type (VvT) ROMK and the 
mutant E104S. The internal solu- 
tion contained 110 mM K + . The 
external solution contained 1 1 mM 
K + plus 99 mM Na + . Lines represent 
fits to Eqs. 2 and 3, with K Mg (0) = 
3.9 mM and z8 Mg = 0.16 (WT), and 



K„ g (0) = 20 mM and z8 : 



Ms 



0.16 



(E104S). (C) Analysis of block by 
3 mM Mg 2 * of WT ROMK and the 
quadruple mutant E92A/ D97V/ 
E104S/E132Q. The internal solu- 
tion contained 110 mM K + . The ex- 
ternal solution contained 1 1 mM K + 
and no Na + . Lines represent fits to 
Eqs. 2 and 3, with KM g (0) = 5.2 mM 
and z8 Mg = 0.20 (WT), and KM g (0) = 
9.3 mM and z8 Mg = 0.17 (E92A/ 
D97V/ E104S/E132OJ. (D) Cur- 
rent traces for WT Kirl . lb in excised 
inside-out patches with different K + 
concentrations on both sides of the 
membrane and 0.5 mM EDTA in 
the pipette solution. (E) Conduc- 
tance as a function of concenuation 
for WT and E92A/ D97V/ E104S/ 
E132Qin the presence and absence 
of 0.5 mM EDTA. Lines represent 
best fits to Eq. 1, with g max = 50 pS ± 1 
and = 11.0 ± 1.3 mM (WT); 
gmax = 53 ± 1 pS and Km = 2.2 ± 0.4 mM 
(WT + EDTA); g max = 47 ± 1 pS and 
K™ = 12 ± 1 mM (E92A/ D97V/ 
E104S/E132Q); and g max = 47 ± 2 pS 
and K,,, = 12 ± 2 mM (E92A/ D97V/ 
E104S/E132Q+EDTA). 
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(Lu and MacKinnon, 1994a). These results contrast 
with the concentration dependence seen in other K + 
channel types. Maxi-K channels have apparent K m 
values of 50-140 mM (Latorre and Miller, 1983). For 
Shaker channels, the value is even higher, around 
300 mM (Heginbotham and MacKinnon, 1993). KcsA 
channels show a more complex biphasic concentra- 
tion dependence with a linear increase in conductance 
at high concentrations of K + (Morais-Cabral et al., 200 1 ) . 

In evaluating the K + concentration dependence, we 
did not substitute other ions for K + . If fixed negative 
charges attract K + to the outer mouth of the channels 
(D'Avanzo et al., 2005; Chang et al., 2010), the change 
in ionic strength would alter the screening of these 
charges, making the measured Km values artificially low. 
On the other hand, if there is a relatively nonspecific 
binding site for cations in the outer mouth, the pres- 
ence of other ions that might compete for this site 
would make the measured Km values artificially high. 
We believe, based on our measurements and as dis- 
cussed below, that the latter artifact is more important, 
so we used a protocol without ion substitution. 

We focused our attention on concentrations between 
5 and 200 mM, where the channels behave as a pseudo- 
single-ion system: 

K c + C < — > K • C < — > Kj + C, 

where C is the channel, and Kj, K • C, and K ; refer to the 
K + ions in the outer solution, bound to the channel and 
in the inner solution, respectively. In such as system, ion 
movement through the channel saturates with an ap- 
parent K m value of k d /k a that reflects the (voltage- 
dependent) rates of association to (k a ) and dissociation 
from (k d ) the binding site. Maximal inward currents are 
determined by the rate of dissociation to the inner solu- 
tion. Under bi-ionic conditions, the relative permeabil- 
ity will reflect channel occupancies by the two species; 
these occupancies will be inversely related to K m . Such 
a model is clearly an oversimplification, as both physi- 
ological and structural measurements indicate that 
the channels are multi-ion pores. However, as shown 



recendy (Nelson, 2011), permeation through multi-ion 
pores can under some conditions still be considered a 
two-step association-dissociation process that follows 
Michaelis-Menton kinetics. 

In any case, the scheme provides some simple approxi- 
mations. From the single-channel currents of 5 pA at 
— 100 mV, k d must be at least 3 x 10 7 s _1 at this voltage. 
If K„ = 10 mM, the value of k a would be >3 x 10 9 M~ ! • s _1 . 
This is close to the maximum rate of ion capture 
by diffusion for a pore with an outer mouth of ^10-A 
radius (Hille, 2001). Cation capture rates by the chan- 
nel may be accelerated by electrostatic attraction to 
fixed negative charges in the channel mouth. 

Rb + and NH 4 + currents 

The apparent Km value for Rb + conduction was 24 mM, 
significandy larger than that for K + . Rb + also has a lower 
maximal conductance than does K + (Fig. 2) , and the 
permeability ratio Pm/Pr based on bi-ionic reversal 
potentials is 0.4—0.6 (Chepilko et al., 1995; Choe et al., 
2000) . This indicates that channel occupancy by K + is 
favored over Rb + . As channel occupancy would in turn 
be inversely related to Kn, as discussed above, the higher 
Km for Rb + is in agreement with the lower permeability. 
The lower maximal currents indicate a lower rate of dis- 
sociation from the channel. 

NH 4 + is a more complicated case. It has a maximal 
inward conductance similar to that of K + , but the per- 
meability ratio is ^0.1 (Chepilko et al., 1995; Choe 
et al., 2000). This could be explained if permeation bar- 
riers for NH 4 + and K + were similar, but that NH 4 + inter- 
acted more weakly than K + with binding sites within 
the channel. Our naive prediction was that the appar- 
ent affinity of the pore for NH 4 + would be lower than 
that for K + . Consistent with this notion, the measured 
K m value for NH 4 + , 14 mM, is indeed higher than that 
for K + . However, the ratio of Km values (Km(NH 4 )/K m (K) 
of <2) is much smaller than that for permeability 
(P k /P N h 4 of ^10). Outward NH 4 + currents apparently 
saturated at lower concentrations. Although these obser- 
vations are difficult to explain with a simple permeation 
model, they can be accommodated by a more realistic 
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Figure 9. Fixed charges in the trans- 
membrane cavity increase affinity for 
K + permeation. (A) i-V relationships 
for Kirl.lb N152D in inside-out patches 
with different symmetrical [K + ]. (B) Con- 
ductance as a function of concentration 
for N152D. Lines represent the best 
fits to Eq. 1, with g max = 39 ± 2 pS and 
K m = 4.7± 1.4 mM. 
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(and more complex) formalism, as shown in Fig. 10 and 
discussed below. 

Outer binding site for cations 

The outer mouth of the pore represents one possible 
location of saturable cation binding. Occupancy of a 
site in this part of the channel could be rate limiting 
for conduction, particularly for inward currents. To 
explore this possibility, we first examined the abilities 
of impermeant ions to block the channel from the 
outside. Small divalent cations (Mg 2+ , Ca 2+ ) blocked 
Kirl.l channels with a small but significant voltage 
dependence (z8 of M).2), consistent with a blocking 
site within the transmembrane electric field. The affin- 
ity of these divalents for the blocking site, assessed as 
apparent K; value at constant voltage, was higher for 
Mg 2+ than for Ca 2+ . A similar block of Kir2.1 channels 
was reported previously (Murata et al., 2002). The im- 
permeant monovalent cations Na + and Li + also blocked 
the channels with lower apparent affinity but with a 
voltage dependence similar to that of Mg 2+ or Ca 2+ . 
Although this could indicate that the monovalent ions 
penetrate deeper into the electric field, it is also consis- 
tent with the idea that monovalents and divalents block 



at the same site, with the voltage dependence attrib- 
utable mainly to the movement of K + ions within the 
field consequent to block. Such a scenario was sug- 
gested for external TEA + interactions of Shaker chan- 
nels (Thompson and Begenisich, 2003). 

Competition between monovalents and divalents in- 
dicates that these ions occupy the same outer site, or at 
least that their binding is mutually exclusive. Based on 
the ability to relieve Mg 2+ block, the monovalent cation 
selectivity for this site is K + > Na + ~ Li + > Rb + ~ NH 4 + . In 
particular, the estimated for K + displacement of Mg 2+ 
was 14 mM, a value not too different from the apparent 
K, n for conduction. In addition, 3 mM Mg 2+ added to the 
outer solution increased the apparent K, n for K + con- 
ductance substantially. These findings are consistent with 
the idea that saturation of K + conductance involves in- 
teraction with this outer binding site. The ability of the 
other permeant ions, Rb + and NH 4 + , to displace Mg 2+ 
was less than that of K + , in qualitative agreement with 
their lower permeabilities and higher K m . 

The SO site, lyingjust outside the selectivity filter of K + 
channels, is a likely structural locus for this outer bind- 
ing site. This site has been identified in crystal structures 
of KcsA channels (Zhou et al., 2001) as well as through 




Figure 10. Kinetic model for permeation through Kirl.l channels. (A) Kinetic scheme with six ion-binding sites and five states. (B) Fits 
of i-V relationships with symmetrical [K + ]. (C) Fits of i-V relationships with symmetrical [NH<*]. (D) Fits of i-V relationships with fixed 
internal [K + ] and variable external [NH 4 + ] . (E) Saturation of K + and NH 4 + currents in the simulated channel. (F) Energy profile for K + 
and NH 4 + movement through the pore of Kirl.l. The major differences are in the well depths of SI and S2. 
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computational modeling (Berneche and Roux, 2001). 
Ions associated with the channel at this location are 
probably partially dehydrated (Berneche and Roux, 

2001) . Occupancy and dehydration may be stabilized 
by electrostatic interactions with negatively charged 
amino acid side chains (Zhou et al., 2001). If this iden- 
tification is correct, this site exhibits a significant though 
modest selectivity for K + over Na + . Another site, termed 
Sexti was also identified as a locus of ion accumulation 
in KcsA channels external to SO (Berneche and Roux, 
2001; Zhou et al, 2001). Blocking at this site is less likely 
to show an appreciable voltage dependence. The block- 
ers are unlikely to bind deeper within the pore, as this 
would predict a stronger voltage dependence, at least 
for divalent cations. If the SI site senses 25% of the 
electric field, block at that site by external Mg 2+ would 
have a z8 value of at least 0.5, and probably more be- 
cause K + would be displaced toward the cytoplasm dur- 
ing occupancy of the site by a blocker. The measured 
was M3.2. On the other hand, the binding of permeant 
cations to the SI site selectivity filter could contribute 
to the ability of monovalent cations to inhibit Mg 2+ 
block at SO, especially in the case of K + . 

The apparent Km for inward currents was weakly voltage 
dependent, at least for K + and NH 4 + , with a higher K m 
estimated when the membrane was hyperpolarized. This 
could arise if the forward rate of exit from the outer 
binding site had a higher voltage dependence than the 
rate of entry onto the site. 

We mutated anionic amino acids in the outer mouth 
that might increase Mg 2+ and/or K + affinities for the 
channel through favorable electrostatic interactions. 
As reported previously (Sackin et al., 2011), the elimi- 
nation of four of these charges (or 16 negative charges 
per tetramer) still resulted in a functional channel with 
conductance and kinetics similar to those of the wild-type 
channel. Mutation of the negatively charged side chain 
closest to the site, E104S, did indeed decrease the affinity 
of the channel for Mg 2+ block; the apparent K(0) in- 
creased from 4 to 20 mM. The elimination of a negative 
charge at an adjacent residue in Kir2.1 had a similar 
effect on Mg 2+ affinity in Kir2.1 channels (Murata et al., 

2002) . The removal of the other three negative charges 
did not further decrease Mg 2+ affinity of the Kir 1.1. The 
apparent K m for K + conductance did not change sig- 
nificantly either with the single E104S mutation or the 
quadruple mutation. In part, this reflects the lower va- 
lence of K + , relative to Mg 2+ , making binding less sensi- 
tive to the electrostatic potential. In addition, however, 
elimination of the screening of the negative charges by 
trace divalent cations by chelation with EDTA revealed 
their impact on K m . 

Inner cation-binding site 

We also examined the effects of adding negative charges 
just inside the selectivity filter within the transmembrane 



cavity. This N152D mutation mimics the natural anionic 
side chain in Kir2.1 that contributes to the strong block 
of this channel by intracellular Mg 2+ and polyamines 
(Lu and MacKinnon, 1994b; Taglialatela et al., 1995; Yang 
et al., 1995). The apparent Km for conductance through 
this mutant was about half that of the wild-type channel. 
Thus, this site can also contribute to the saturation be- 
havior of inward rectifiers. 

A kinetic model 

Our simulations addressed the question of whether 
standard models of permeation through K + channels 
with designated ion-binding sites identified by x-ray 
crystallography could explain the observed behaviors. 
We could identify a set of kinetic parameters that ac- 
counted for many of the channel properties, including 
the low Km values for both K + and NH 4 + , the high K + /NH 4 + 
permeability ratio despite similar maximal conductances, 
and the increase in inward NH 4 + conductance with K + 
on the opposite side of the membrane. The model is 
based on discrete-state kinetics and contains a minimal 
number of sites consistent with structural information. 
We do not claim that it is unique in describing the ex- 
perimental results. However, we believe that the sim- 
ulations show that our data can be explained without 
postulating any new or unusual features of the channel 
or a qualitatively different interaction of NH 4 + and K + 
ions with the channel. 

Analysis of this model helps to explain a major con- 
clusion from experimental findings: two observable 
parameters related to the strength of interaction of 
permeant ions with the pore — bi-ionic permeability and 
apparent affinity — can be dissociated. The relative per- 
meability of ions depends on the energy of binding to all 
the sites within the multi-ion pore. Saturation of inward 
conduction through the channels reflects the occupancy 
of a cation-binding site outside the selectivity filter and is 
not necessarily correlated with permeability. The similar 
Km values for permeant cations reflect the relatively non- 
specific cation binding to this part of the channel. 
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